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ABSTRACT 

The electrical resistivity of the palladium-deuterium system was measured to a 
deuterium- to palladium -atom ratio of 0. 9 at temperatures of 273, 77, and 4.2 K. The 
resistivity ratio p(x)/p(0) was plotted as a function of the atom ratio x at 273 and 
4.2 K. A modification of Mott’s model for the resistivity of transition- metal alloys was 
used to calculate the structural resistivity. A good fit to the data at 4.2 K was obtained 
by assuming that the number of d-holes per palladium atom takes on the value of 0. 55 to 
0.60 on addition of deuterium. 
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ELECTRICAL RESISTIVITY AS A FUNCTION OF DEUTERIUM 
CONCENTRATION IN PALLADIUM 
by Gust Bambakidis, Robert J. Smith, and Dumas A. Otterson 
Lewis Research Center 

SUMMARY 

The electrical resistivity of the palladium -deuterium (Pd-D) system was measured 
to a deuterium- to palladium -atom ratio of 0. 9 at temperatures of 273, 77, and 4. 2 K. 
The resistivity ratio p(x)/p(0) was plotted as a function of the atom ratio x at 273 and 
4. 2 K. A modification of Mott's model for the resistivity of transition-metal alloys was 
used to calculate the structural resistivity. A good fit to the data at 4. 2 K was obtained 
by assuming that the number of d-holes per Pd atom takes on the value of 0. 55 to 0. 60 
on addition of D. 


INTRODUCTION 

At room temperature, deuterium (or hydrogen (H)) occupies the octahedral interstitial 
sites of the face-centered-cubic Pd lattice (refs. 1 and 2). The octahedral sites are oc- 
cupied for both a-Pd and j3-Pd, with the D/Pd and H/Pd atom ratios ranging from 0 to 
0. 72 (where a-Pd and /3-Pd have lattice constants of 3.86xl0”^and 4. 02x10”^ cm, respec- 
tively). Denoting the D/Pd or H/Pd atom ratios by x reveals that only a -Pd exists for 
0 < x < 0. 02. For 0. 02 < x < 0. 60, a-Pd and /3-Pd coexist, and only /3-Pd exists for 
x > 0. 60. If p(x) is the resistivity of PdH (or PdD ) and p(0) the resistivity of annealed 
Pd, the resistivity ratio [p(x)/p(0)]2gg K and its derivative increase for increasing x 
at 0. 02 < x < 0. 70. The resistivity ratio is a maximum for x ~ 0. 76, and the ratio 
along with its derivative decreases for 0. 76 < x ^ 0. 88 (ref. 3). In the PdD x system, 
the available data show that the ratio Ip(x)/p( 0 )] 2 98K and its derivative increase from 
x 0. 02 to x = 0. 67 with no higher values of x being shown (ref. 4). For pure a-Pd 
(x < 0. 02), the derivative of the resistivity ratio [p(x)/p(0)] 2 gg K is larger than that for 
the a-Pd + /3-Pd phase for either PdH x or PdD x . 



The objective of the present work is to extend the resistivity data for PdD x beyond 
x = 0. 67 to see if, in the plot of p(x)/p(0) as a function of x, a maximum occurs similar 
to that found for PdH x at 273 K. The thermal contribution to the resistivity was mini- 
mized by obtaining data at 4. 2 K, and then a modification of Mott's model for the resis- 
tivity of transition- metal alloys (ref. 5, p. 296) was used for the analysis of the struc- 
tural contribution. The resisitivity results and analysis presented herein for the PdD 
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system should be qualitatively applicable to the PdH system since these systems are 
quite similar in other resistive and crystallographic properties thus far investigated 
(refs. 1 to 4 and 6 to 8). 


EXPERIMENTAL PROCEDURE 

The experimental data were obtained at 4. 2, 77, and 273 K. Obtaining the 4. 2 K 
data necessitated rapidly cooling or quenching the PdD x wires from 273 to 4.2 K to 
prevent the migration of deuterium ions from their octahedral sites (0, 0, 1/2; 0, 1/2, 0; 
1/2, 0, 0; . . .) to their tetrahedral sites (1/4, 1/4, 1/4; 3/4, 3/4, 3/4; . . .). For 
PdH x , an analysis of the data for neutron diffraction and that for resistivity as a function 
of temperature p(T) showed that up to 25 percent of the hydrogen may migrate to 
tetrahedral sites. These octahedral-tetrahedral transitions occur for PdH at 

X 

0. 45 < x < 0. 75 if the PdH wires are cooled without quenching (refs. 6 and 7). In the 
PdH x system, the octahedral-tetrahedral transitions may begin well above 100 K and 
continue at least to the neighborhood of 50 K (at which point a maximum in the resistivity 
and specific heat occurs) as the temperature is lowered from 273 to 4.2 K (ref. 6). 
(Whether or not additional transitions occur below these maximums is not known.) Al- 
though no neutron diffraction data are available for PdD , the resistivity data of PdD 
are almost identical to those obtained for PdH between 273 and 4. 2 K (ref. 8). How- 

X 

ever, at a given temperature, the octahedral-tetrahedral transitions take much longer in 
the PdD system (up to several hours) than in the PdH system (not more than a few 
seconds). From these data for PdD , one can conclude that rapid cooling of the PdD 
specimens from 298 to 4.2 K, as employed in the present experiments, was sufficient to 
prevent a significant amount of deuterium migration to the tetrahedral sites. 

Wire specimens of 99.995 atomic percent Pd (manufacturer's stated purity) and 
0. 254-millimeter diameter were annealed by joule heating in a vacuum of 5x10 torr at 
approximately 1200 K for at least 0. 5 hour. Measurements of initial resistivity p(0) 
were made at 4. 2, 77, 273 K by using the conventional four-probe potentiometric tech- 
nique in which a phenol formaldehyde specimen holder utilized phosphorbronze spring 
contacts for the current and potential contacts to the specimen. Six potential contacts, 
spaced 1 centimeter apart between the current contacts, permitted five consecutive 
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measurements along the wire. This spacing allowed the detection of resistance irregu- 
larities due to preabsorption conditions and/or nonuniform absorption along the wire 
after deuterium absorption. The measuring currents used were 10 milliamperes in all 
cases 

Deuterium absorption was accomplished electrolytically with a graphite anode and a 
Pd-wire cathode in a 1-normal solution of deuterium sulfate (DgSO^). The current was 
5 milliamperes. The anode and cathode were separated by a porous alumina crucible, 
as shown in figure 1 . The separation was to prevent oxidizing agents that might be formed 
at the anode from reaching the cathode, which would permit the Pd to absorb higher con- 
centrations of deuterium. Immediately after the absorption process, the wires were 
returned to the phenol formaldehyde holder for resistivity measurements at 4.2, 77, 
and 273 K to obtain p(x) or the resistance R(x). 

For determination of the D/Pd atom ratios, each specimen was sectioned into 
five 1- centimeter pieces that corresponded to the positions on the phenol formaldehyde 
specimen holder. Then the deuterium was removed by heating each 1-centimeter section 
separately to 573 K in an evacuated system (~10”° torr). Tests showed that less than 
1 percent of the total gas concentration generally remained in the specimens. The spec- 
imens were weighed to three significant figures. The deuterium was expanded into and 
analyzed with a mass spectrometer, which is capable of 1 percent precision. The D/Pd 
atom ratios were then compared with their respective resistivity measurements. 
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Figure 1. - Container used for deuterium absorption by palladium. 
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RESULTS 


The results of resistivity measurements made at 4.2 and 273 K are shown in fig- 
ure 2. As seen in figure 2(a), [o(x)/p(0)] 27 g K and its derivative increase with x for 
0. 02 < x < 0. 65, and [p(x)/p(0)] 217 g K reaches a maximum in the neighborhood of 
x ~ 0. 72 and then decreases as x — 0. 90. At 4. 2 K, the maximum in [p(x)/p(0)]^ 2 k 
appears to be near a value of 0. 68 for the D/Pd atom ratio. 



(a) Experimental curve of resistivity ratio at 273 K as func- 
tion of deuterium- to palladium-atom ratio. 



(b) Resistivity ratio at 4. 2 K. Curves calculated for num- 
ber of d-band holes n h equal to 0.36 and 0.55. 

Figure 2. -Resistivity ratio as function of deuterium 
concentration in palladium. 
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Some of the scatter in the data may be attributed to two sources: (1) some hydrogen 
contamination, which varied between 3 and 8 percent of the total gas concentration, and 
(2) some evolution of the absorbed deuterium, for x > 0.7> between the resistivity mea- 
surements and mass spectrographic analysis. The time period for which evolution could 
take place was less than 20 minutes. (There is no evidence that evolution occurred at or 
below 77 K. ) The possible gas loss for this time period is approximately a D/Pd atom 
ratio of 0.025, which is equivalent to a 2-percent change in p(x)/p( 0) at 273 K. The 
data shown in figure 2 include both the deuterium and the hydrogen impurity in x 
(i.e. , (D + H)/Pd). 


DISCUSSION 

In figure 2(a), an experimental curve was drawn through the data points, whereas the 
curves in figure 2(b) were calculated for the values 0. 36 and 0. 55 for the number of 
holes n^ per Pd atom in the d-band. (The results for 77 K are not given because some 
octahedral-tetrahedral transitions occur during measurement. ) The difference in the 
data at the two temperatures shown is attributed to the thermal contribution to the resis- 
tivity. At 4. 2 K, both the thermal resistivity and that arising from lattice defects are 
expected to be small, so that their variation with deuterium concentration can be 
neglected. If x denotes the D/Pd atom ratio, the resistivity ratio p(x)/p( 0) at this 
temperature can therefore be written 


£« = ( 1 ) 

p( 0) P( 0) 

where pjj(x) is the contribution from disorder scattering by the deuterium randomly 
distributed among the octahedral sites in the host face- centered- cubic Pd lattice. This 
quantity is also the resistance ratio because the change in sample dimensions due to 
lattice expansion on absorption of D is negligible here. The average value of p(0) for 
the samples used was determined to be 0. 105 microhm- centimeter. 

Magnetic susceptibility (ref. 5, p. 189) and electromigration (ref. 9) studies on the 
Pd-H system indicate that H exists in the Pd lattice as strongly screened positive ions, 
and presumably the same situation exists for D. The electronic structure of pure Pd 
was described by use of the Mott model (ref. 5, p. 189), in which the band structure is 
assumed to consist of a broad free-electron-like s-band overlapping with a narrow 
d-band, so that the Fermi level occurs near the top of the d-band. Thus, there are 
electrons in the s-band and an equal number of holes in the d-band. The assumption is 
made that initially the D atoms contribute their electrons to fill the holes in the d-band. 
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In a covalent bond picture, the electron is thought of as being shared with neighboring Pd 
atoms, which results in a strong screening of the D ion. The D ions, therefore, scatter 
a conduction electron only weakly, giving rise to a small initial slope (the change in slope 
on appearance of the /3-phase, observed at higher temperatures, is greatly diminished at 
4.2 K, which suggests that its origin is thermal rather than structural). As the holes in 
the d-band are filled, the electrons are instead contributed to the s-band, whose states 
are delocalized and, hence, do not screen as strongly. This lack of screening results in 
the sharp rise in resistance. Thus, the resistance is viewed as arising from scattering 
by three types of local structure, denoted by I, n, and III in the potential energy diagram 
of figure 3. Each of the three types is distributed randomly throughout the system. The 
screened coulomb potentials about the D ions have screening lengths a^ and a g , with 

a d < v 



Figure 3. - Potentials associated with the three types of local structure 
assumed present in palladium-deuterium system. I, palladium atom asso- 
ciated with strongly screened deuterium ion, screening length a d ; II, 
palladium atom associated with weakly screened ion, screening length a s ; 
III, palladium atom associated with a vacancy. 


If the current carried by the holes is neglected, the use of the Boltzmann equation 
and the assumption of an isotropic relaxation time r lead to the usual expression 


P = 



N e e 2 r 


( 2 ) 


where m* is the effective mass and N the number of conduction electrons per unit 
c e 

volume. The inverse of the relaxation time can be written 
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( 3 ) 
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S-S 
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s-d 


However, the contribution from s-d scattering, proportional to the density of states at the 
Fermi level, goes to zero as the d-band fills. We therefore feel that s-d scattering 
cannot play a significant role in the sharp rise in p(x)/p(0) with x and neglect it 
entirely. 

When 1/t is evaluated in Born approximation (see appendix), the result for 

S-S 

P D (x) is 


2 m e 2e2 

P D ( X ) = 

D 3 * h 3 




_ f K 2a d) x d( 1 - x d) 


+ f (4k 2 a 2 )x s (l - x s ) + 2g(k t ,a d ,a s )x d x s J 


(4) 


The quantity fi Q (x) is the volume per Pd atom which is weakly dependent on x because 
of lattice expansion. Because the exact variation of lattice parameter with concentration 

O 

is not important here, a simple linear expansion of the lattice parameter from 3. 88x10 

O 

centimeter at x = 0, to 4.02x10” centimeter at x = 0. 60, is assumed. The quantity 
n g (x) is the number of electrons per Pd atom in the s-band. The Fermi radius is 
given by 


kf = 


37Tn e (x) 


n o (x) 


1/3 


while f and g are defined by 


and 


f (4kj a* 2 ) = In (l + 4k^a 2 ) - 


4l£a 2 


1 + 4k 2 a 2 


g(kf,a d , a s ) - 


4 !n (l + 4lja 2 ) - a 2 In (l + 4k f 2 a^ 

2 2 
a - a, 
s d 
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The coefficients x^ and x g denote the number of D atoms per Pd atom contributing an 
electron to the d-band and s-band, respectively. Both depend implicitly on the total 
concentration x. 

Expression (4) cannot be expected to be valid in the upper range of concentration (at 
and above the maximum in [p(x)/p(0)] 4 b ecause the development of long-range order 
near x = 1, where the system takes on the sodium chloride structure, was not taken into 
account. The expression, therefore, overestimates the resistance in this range. If at 
x = 1, type I and n structures are regarded as constituting independent super lattices, 
then the resistance for x < 1 will arise from scattering by type II and in structures 
only. Proceeding as before in the evaluation of 1/r results in the high- charge 
expression (see appendix) 



*2 2 
m e e 


IT 


n 0 (*) ■ 



(x - n h )(l - x) 
! -n h 


( 5 ) 


where x > n^. The parameter n^ is the number of holes per Pd atom in the d-band at 
the initiation of charging; hence, it is also the number of octahedral sites per Pd atom 
available to the strongly screened D ions. The assumption of independent superlattices 
means that equation (5) underestimates the resistance in the range x > n^. 

The terms x^, x g , and n g are expected to have the qualitative behavior shown in 
figure 4. Note that x . + x 0 = x and that n (x) = n (0) +x . The concentration x, is 
that at which the d-band fills, and experiments by others (ref. 5, p. 316) on palladium - 
noble-metal alloys indicate x^ to be in the region of 0. 55 to 0. 60. On the basis of the 



Figure 4. - Expected total -concentration dependence of 
quantities x d , x s , and n p entering into equations (4) 
and (5). 
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de Haas - van Alphen measurements of Vuillemin (ref. 10) on pure Pd, n e (0) = 0. 36. If 
this value is also used for the number of holes, the upper curve in figure 2(b) is obtained. 
The dashed line indicates the continuous curve that would result from a more exact 
treatment of the onset of long-range order. The effective mass used, which was assumed 
independent of composition, was an average band mass of 0.49 m g obtained from k^ and 
the Fermi level for pure Pd. The quantity was taken as 0. 462 rydberg by using 
Segall's (ref. 11) band- structure calculations for copper and the similarity in the band 
structures of Pd and Cu (ref. 10). The quantity x^, was taken as 0. 60, and a^ and a g 
were determined by fitting equations (4) and (5) to the initial slope and final (extrapolated) 
slope, respectively. The fit over the entire concentration range is certainly not quan- 
titative, but if the number of holes in the d-band is assumed to increase from 0. 36 on 
charging, a much better fit is obtained, as shown in figure 2(b) by the curve for 
n h = 0. 55. Values for n h in the range 0. 55 to 0. 60 are consistent with magnetic sus- 
ceptibility data obtained by others on Pd-H (ref. 5, p. 316). The values for the screening 
lengths are a^ = 0. 26 atomic units and a g = 1.11 atomic units (at x = 1). This value for 

a_ is somewhat larger than the Thomas-Fermi value of 0. 74 atomic units, 
s 

The authors are aware that the assumption of a spherical s-electron Fermi surface 
neglects the appreciable anisotropy, which extends to about 30 percent (ref. 10). How- 
ever, the Mott model, with the interpretation of the de Haas - van Alphen work that 
n g (0) = n^O) = 0. 36, appears to work well in explaining the transport properties of the 
Pd - noble-metal alloys (ref. 12). The apparent failure of the rigid-band model for Pd-D 
indicated herein suggests that further work on this system would be of interest. In parti- 
cular, a band- structure calculation could be done for the ordered structure x = 1 to 
obtain n g (l). The rigid-band model predicts one s-electron per Pd atom at D/Pd or at 
H/Pd = 1, whereas the present analysis would indicate 0. 81 electron per atom. In con- 
junction, a de Haas - van Alphen experiment could be conducted, which would be depend- 
ent, however, on overcoming present experimental difficulty in consistently obtaining 
high D/Pd ratios. 


CONCLUSIONS 

The direct-current electrical resistivity of the palladium-deuterium (Pd-D) system 
was measured to a deuterium- to palladium-atom ratio of 0. 9. Analysis of the resistivity 
as a function of atom ratio at 4.2 K indicated that on addition of deuterium, the number of 
holes in the d-band (per Pd atom) takes on the value of 0. 55 to 0. 60, compared with the 
value of 0. 36 for pure Pd, which is based on recent experimental and theoretical work by 
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others. This result suggests that the rigid-band model is not valid for the Pd-D inter- 
stitial solid solution, in contrast to the Pd - noble-metal alloy systems. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 16, 1968, 

129-03-15-01-22. 
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APPENDIX - EVALUATION OF INVERSE OF RELAXATION TIME 


The inverse of the relaxation time t for scattering of an s-electron at the Fermi 

S — s 

surface back into an s-state by a single scatterer is (ref. 5, p. 247) 


1 


T 


s-s 



dS(l - cos 



JL 

27m 



d0(l - cos 9) sin 


(Al) 


The integral is taken over all scattering angles 9 between the initial state k and final 
state k'. The matrix element v k , k is 

f dr ^J,(£)v(r)^ k (r) 

where is the initial state in the presence of the scattering potential v(r), and <p k , 
is the final state in the absence of v. The s-band wave functions are assumed to be ~ 
plane waves. In the Born approximation, i^ k is taken equal to cp^. in evaluating v k , k . 
Hence, — 


v k , k =i / dre-^’-B-IvW (A2) 

« •'n ~ 

Since a perfectly periodic potential has no resistance, the contribution from each 
type of local structure of figure 3 is considered to arise from the difference between the 
actual potential and the average potential there. The average potential is 

v av = [ v Pd + v D (a d>]*d + [v pd + v 0 (a s )]x s + v pd (l - x d - *„) 

Hence, the scattering potentials are 
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Type I: 


V (r) = v - v ay 

= [ v Pd + v D< a o3 ' IjPd + X d v D< a d> + x s T D< a si] 
= (1 - x d )v D< a d> - x s v D< a s> 

Type II: 

V (r) = (1 - x s )v D (a s ) - x d v D (a d ) 

Type HI: 


v(r) = -x d v D (a d ) - x g v D (a s ) 

The total contribution to l/r_ _ from each type of scatterer consists of a term of the 

s — s 

form of equation (Al) with the appropriate scattering potential, weighted by the number 
of each type. Denoting matrix elements by (( }) and integration over the scattering angle 
by ( ) gives 

— = ““ k f(~) 1 {J^dK^I 2 x d^ " x d> + 1 < y D ( a s)> I 2 x s (1 “ x s) 

T s-s 27m Vdk/ k=k f 


+ 2 Re[<v D (a d )><v D (a g ))*] x d x g | (A3) 

where Np^ is the number of Pd atoms in the system. The appropriate integrals can be 
evaluated in a straightforward manner, the result being equation (4) of the text. 

In the high- charge region, the resistance is assumed to arise mostly from the pres- 
ence of vacancies in a superlattice of type II structures. The scattering arises from 
fluctuations in an average potential: 


v av = v Pd + 


S v D< a s> 


1 - n. 


The scattering potentials are now 
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Type II: 


v(r) = 





Type HE: 


v(r) = - 


1 - n v 


v D ( a s ) 


There are x type n scatterers and 1 - x type III scatterers per Pd atom. Therefore, 


1 _ ^ N Pd k 2/dE\~ 1 


T S-S ^ 


) 


Kv ^! 2 




+ ! < v D (a s» i 2 



(A4) 


Now, x^ + x g = x, and from figure 4, x^ ~ for x » n^. Thus, using x g = x - n^ 
results in 


s-s 


^my 1 i < v p( a s )) i 2 (x ~ ° h)(1 

2l * Vdk/ k=k l-n h 


x > n. 


(A5) 


Substituting the explicit expression for | (vj-j(a g )> | 2 leads to equation (5) of the text. 
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